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Abstract 

The present paper proposes Hybrid Fractal image 
compressions (HFIC) for grey and Color images. The 
HFIC method is carried out for the variable size range 
blocks. The sample images here are divided into blocks by 
considering minimum and max size ranges. This method 
transforms the RGB color image into three planes and then 
divide them into variable ranges. To achieve fractal 
compression, Relative transformation and entropy coding 
is applied. The image is reconstructed by applying the 
inverse transforms and iterative functions. The results of 
the present analysis demonstrate that, for the hybrid fractal 
compression scheme with variable range method applied 
to color and gray scale images like Lena, medical MRI etc, 
show higher CR and PSNR values compared to fixed 
range block size and other existing methods. The results 
indicate that both color and gray level images with R max = 
16 and R m in = 8 shows higher Compression ratio(CR) and 
good Peak Signal to Noise Ratios (PSNR). 

Keywords: Hybrid fractal image compression (HFIC), 
Minimum range Block, maximum range Block, Relative 
transformation and entropy. 

Nomenclature 

HFIC Hybrid Fractal Image Compression 
MRI Magnetic Resonance Imaging 
FIC Fractal Image Compression 
PSNR Peak Signal to Noise Ratio 
CR Compression Ratio 

MFCS A Modified Fractal Coding on Spiral Architecture 
GA Genetic Algorithm 
D-Block Domain Block 
R-Block Range Block 

1. Introduction 

A fractal can be defined as a geometrical shape that is self- 
similar i.e. it has parts that is similar to the whole. Fractal- 
based image compression (FIC) techniques exploit 
redundancy due to self-similarity properties in images to 
achieve compression [1]. Fractal image compression 
schemes for image compression have been developed and 
implemented. 


Generally, Image compression means minimizing the size 
in bytes of a graphics file without degrading the quality of 
the image to an unacceptable level [1]. The main motto of 
reduction in file size is more images to be stored in a given 
amount of disk. Another advantage is reduces the time 
required to upload and download the images over an 
internet. The main principle of the image compression is 
that the neighboring pixels are correlated and therefore 
contain redundant information. The foremost take of the 
Image Compression is to find the correlated representation 
of the image[2]. 

The idea of Fractal Image Compression (FIC) was 
originally introduced by M.F.Bamsley and S.Demko [3]. 
Fractal image compression that uses the characteristics of 
existing self similarity [4, 5] within images is a suitable 
method for coding an image. Jacquin introduced fractal 
image encoding based on PIFS [6] later. Fractal Image 
Compression (FIC) is the most widely used applications 
such as image retrieval, texture segmentation, feature 
extraction and image signature in the field of image 
processing [7]. Fractal compression allows fast 
decompression but has long encoding times. The most 
time consuming part is the domain blocks searching from 
each range [8, 9]. 

2. Literature Review 

From the past few years, so many fractal image 
compressions techniques are proposed by many researches 
and authors. In fractal image compression, the image is 
separated into various domain blocks with subjective size 
ranging from 4^4 to 16x16, or more. At that point, the 
picture is divided again into range blocks of short size than 
that of the block domain [10]. For each domain range 
match two transformations are required, a geometric 
change which maps the domain to the extent and an 
relative transformation that modifies the intensity values 
in the domain to those in the extent [11]. The fractal 
compression technique as explained in [12], [13], [14] is 
basically a search process consists of partitioning the 
image into sub images and search for parts of the images 
which are self similar. The various partitioning schemes 
are compared in [15]. The algorithm used for encoding is 
the Partitioned Iterated Function System [16] compared 
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with other image compression methods, hybrid fractal 
with binary plane algorithm [17]. 

M.Ismail et al [18] proposed a fractal image compression 
system focused using range size. Rama Krishna et.al [19] 
executed lossless compression using Huffman coding 
working straightforward with quantized coefficients. 
D.Venkatasekhar et.al [20] proposed FIC method for 
image compression based on Polynomial hybrid Wavelet 
and particle Swarm optimization. This method use DCT 
and DWT techniques but it takes most computing time for 
domain blocs. The method is applied on only three 
categories of images such as parrot, car and flower. When 
this method is applied on other categories of images good 
and comparative results are not obtained. 

Wei Hua Rui et al., proposed a progressive texture 
compression framework [21] to reduce the memory and 
bandwidth cost by compressing repeated content within 
and among large-scale remote sensing images. M. Ashok. 
T. Bhaskara Reddy et al proposed, Color Image 
Compression based on Luminance and Chrominance[22] 
Nileshsingh V. Thakur, Dr. O. G. Kakde [23] proposed the 
approach Modified Fractal Coding Algorithm for Grey 
Level Images on Spiral Architecture (MFCS A), 
composing with the one-plane image using the pixefs 
trichromatic coefficients. One-plane image in traditional 
square structure is represented in Spiral Architecture for 
compression. On this Spiral Architecture image, proposed 
modified Fractal grey level image coding algorithm 
(MFCS A) is applied to get encoded image. 

The research has been carried out to achieve the fast speed 
of fractal encoding and also high quality of reconstructed 
image [24, 25, 26, 27, 28]. From the above literature 
observed that no single method is applied for all categories 
of images. 

A.R. Nadira Banu Kamal et.al [29] proposed an efficient 
search of the domain pools for color image compression 
using Genetic Algorithm (GA). The proposed method 
reduces the coding process time and intensive computation 
tasks. Image quality factor, compression ratio and coding 
time are analyzed in this approach and comparative results 
are obtained but this method is also suitable for only 
limited categories like facial, building and vegetable 
images. K. Raja Kumari et.al [30] proposed a fractal image 
compression technique based discrete wavelet transform is 
followed by Huffman Run length Encoding. This method 
is mainly focused on improving the image quality after 
compression. This method is also not suitable for all 
categories.. Bakul Pandhre et.al; [31] proposed an 
approach to reduce the encoding time of an image by 
classifying the blocks according to an approximation error 
measure. This method is also applicable not applicable for 
many categories of images. 

From the above literature survey we observed that no such 
method is available to compress the all categories of 
images by using FIC technique with good image quality 
after compression. This problem is addressed here and a 
method is developed called Hybrid Fractal image 
compression (HFIC) to compress major categories of 
images like standard pepper , lena medical MRI , Internet 
images etc. 

The paper is organized into the following sections. Section 
3 describes the preliminaries of basic fractal image 
encoding and decoding approach and section 4 emphasizes 


on the Hybrid Fractal image compression. Results and 
discussion are focused in section 5 and conclusions are 
given in section 6. 

3. Fractal image encoding and decoding 
process 

The essential steps of fractal image compression (FIC) are 
as follows 

Step 1 : Partition original image into non-overlapped little 
sub windows. Each sub window is called Range-R block 
Step 2: For each R-block, discover an overlapping image 
block i.e. Domain block or basically called D block. The 
D block which is most similar to current R-block under a 
certain change 

For example, if the size of the image is 256x256 which is 
to be encoded having 256 gray shades. If the size of the R- 
block is considered as 8x8 then the available R-blocks in 
a entire image is 1024. These entire R-blocks are 
composite into R-pool. Suppose D-block is four times 
larger as R-block viewed then the quantity of accessible 
D-blocks is (256-2*8+l) 2 =58081 and all these D-blocks 
are composite D-pool. 

Step 3: For every R-block, locate suitable D-block from 
the D-pool which is most parallel to it. 

The concrete steps for discovering the suitable D-block 
from D-pool which is most parallel to corresponding R- 
block are 

(i) Crop the D-block size to the size of R-block by Shrink 
operation and, marked the D-block as D’ -block using Four 
neighborhood regional method shrink operation. 

(ii) Swap the D’ -block by using eight relative 
transformations matrices which proposed by Jacquin. 
Shown in figure 1 . 


!) 

MA o) 

Ti -T o 1 -°i) 

*-(? -;) 

M'o 1 ?) 

*-G -°i) 

M-°i “o 1 ) 

j) 


Figure 1. Relative Transformations Matrices 

(iii) Eight relative transformation blocks are generated for 
each D’ -block and these blocks now composite new D- 
pool. 

(iv) Match up to each R-block with the blocks which are in 
D’-pool and obtain nearly all similar blocks. The matching 
can be measured with normal clash MSE. Generally the R- 
block and D-block are vectors thus, the MSE can be 
measured by using the following equation 1 . 

MSE =//{X k *D k +Y k *1-R k )// (1) 


Where X and Y are the coefficients and should have the 
following values for making the MSE values minimum. 
The coefficients are considering form well known least 
square method. The Xk and Yk are represented in the 
equations 2 and 3 respectively. 


Xk = 


B 2 <D k> R k > ~<D k ,IXR k ,I> 
B 2 <D k ,D k > -<D k ,I > 2 


( 2 ) 
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y _ <R k J> <D k ,l> 

Y k - ( 3 ) 

The meaning of the above step is choosing a block with 
the smallest approximation of D-block’ s eight 
transformations matrices for R-block with least square 
approximation. Then compression relative transformation 
W. W = {Dk(x,y), Tk, Xk, Yk). of each R-block is expressed 
in the equation 4 Where Dk(x,y) represents the D’-block’s 
staring position and Tk represents the corresponding 
relative transformation matrix. 

(v)Store up the compression relative transformation W for 
each R-block. All clamping relative changes constitute the 
entire image's fractal code. 

The figure 2 represents the schematic diagram of the 
conventional fractal encoding process when applied on 
Lena image. 



choose the best 
approximate result 

For each D block, use least square 
method to approximate R block 


eight 

transformed 

block 


Figure 2: Schematic diagram of the fractal encoding process 

3.1 Fractal image decoding process 

The fractal decoding procedure is quite simple, 
that is, for all R blocks, act the corresponding relative 
transformation on any one initial image, iterate for several 
times(commonly ten times) , according to fixed point 
theorem, this iteration process would converge to a fixed 
attractor, and this attractor is the decoding image. 
Adaptive threshold quad tree fractal compression 
approach is the same as the basic fractal image encoding 
method, the main difference is, when doing image match, 
we don’t find the global best-match block in D-pool, and 
instead, we use an adaptive threshold to judge whether the 
current R-block matches D-block, if the two’s match error 
MSE is smaller than threshold, we can consider that 
current D-block matches R-block, otherwise, the two are 
not matching. The threshold is proportional to the image 
block’s variance. When block’s variance is larger, i.e., 
block has a high complexity, and has much high-frequency 
signal, we can allow a greater distortion in recovery image 
block. When block’s variance is smaller, i.e., block has a 
low complexity, and has much low-frequency signal, we 
must limit the distortion to a certain level. 

Points observed from FIC approach 

(i) Fractal image compression approaches do not consider 

the pi^ure’s content. 


(ii) FIC divides image into regular square image sub block. 
This effects visually recovered image when the main 
image block is larger in size, reducing SNR. But when 
image block has a smaller size, it can increase D-blocks’ 
count, thus increasing the encoding time and reducing 
compression ratio. Inspired by this, it is proposed to use 
some statistic information of image block (such as gray 
average, variance, fractal dimension) to speed the 
matching process reducing time complexity. 

From the above discussion, it is clear that in 
fractal image compression, the encoding time and 
intricacy is high furthermore distinctive techniques are 
utilized for different images. To overcome these issues, an 
efficient hybrid scheme focused around variable range 
block with quadrant methodology is proposed for 
encoding. 

4. Proposed Hybrid Approach 

The main objective of proposed work is to increase the 
effectiveness of encoding and decoding process in FIC and 
also compress the RGB images instead of grey level 
images. The proposed hybrid algorithm enhances the 
fractal image coding and attempts to improve and increase 
the compression ratio while maintaining a high 
reconstructed image quality. 

FIC is basically dividing the original image into 
non-overlapping regions called R blocks and overlapping 
regions called D blocks. For each R block, the greatest 
matching D block must be found by relative changes Wi is 
of the form as follows in equation (4). 



-x- 


"CZj 

bi 

0- 


-x- 


[Gl 

Wi 

y 

= 

G 

d t 

0 


y 

+ 

ft 


z 


.0 

0 

Si. 


z 


Oi. 


In above relative changes w if Si manages the contrast and 
Oi controls the brightness and a u b i} c iy d h e if f means the 
eight symmetries like for example, identity, rotation 
through +90°, rotation through +180°, rotation through - 
90°, reflection about mid-vertical axis, reflection about 
mid-horizontal axis, reflection about first diagonal and 
reflection about second diagonal. Figure 3. shows the 
proposed hybrid fractal image compression system for 
color images and proposed algorithm is shown below 


Partitioning 



Figure 3: The Schematic diagram of the Proposed Method 

4.1 The Hybrid Fractal Encoding and Decoding 
Algorithm: 

The algorithm steps are as follows 
(i)Separate R, G, and B components of an image. The 
image is converted to three 2-D matrices i.e. Red plane 
image, Green plane image and Blue plane image. 
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(ii) Each plane image will be partitioned into non- 
overlapping blocks. These are called Range block or 
simply R-block. Initially select the R-block size ranges. 
The maximum Range block of size (R ma x) is 16 or 8 and 
minimum Range block of size (R m in) is 4 or 8. R-blocks 
are compared with domains from the domain pool, which 
are twice the range size. 

(iii) The domain block with the size of window KxK are 
sliding over the entire 3 separated images in steps of K/2 
or K/4 known as lattice. The pixels in the domain are 
averaged in groups so that the domain is reduced to the 
size of the range and applying relative transformation. 

(iv) After partitioning and transformation, the fractal 
encoding process is the search of suitable candidate from 
all available blocks to encode any particular range block. 

(v) Attempts to improve encoding speed involves 
classification of sub-image into upper right, upper left, 
lower right and lower left quadrants shown in Figure. 4. On 
each quadrant compute values proportional to the average 
intensities of the 3 plane images. They will follow one of 
the three ways as canonical ordering . 



Alpha N U 

Figure 4: Classification Scheme of the Sub Image by Canonical 
Ordering 

They are (a) Major Class 1 : vi>V 4 >V 2 >V 3 

(b) Major Class 2: vi>V 3 >V 2 >V 4 

(c) Major Class 3 :vi>V3>V4>V2 

(vi) Addition to the above three major classes, there are 24 
different subclasses for each major class. So totally 
domain and range blocks are represented in 72 classes. In 
coding process any range block is mapped to the domain 
blocks and using of the entropy coding to achieve fractal 
compression for each plane image 

(vii) Record the fractal decoder to re-construct the image 
and calculate CR and PSNR. 

5. Results and Discussions 

To find the effectiveness of the proposed method 
illustrated by means of the experimental results. The 
proposed approach was executed in Matlab 13b on Intel 
core i3 processor with 2 GB RAM and the proposed 
method was evaluated on both gray and color images. The 
Proposed technique is applied on a set of color images of 
different categories with the standard size512><512. The 
test images used in the experiments are the color Standard 
Lena image of size 512x512, Medical images with the 
resolution of 512x512 is obtained from Hospital, Pepper 
image with size of 512x512, Zelda image from goggle 
with the size of 512x512 and God hills images form 
Google with the size of 600x800. In our experiment all 
images consider the standard size 512x512. The Quality 
of the re-constructed images was determined by measuring 
the Peak Signal to Noise Ratio (PSNR) value and the 
Compression Ratio(CR) for both color and grey level 
images. 

By using the variable range block size for three cases 
namely (a) Rmax - 1 6 and Rmin — 4 (b) Rmax — 1 6 and Rmin — 
8 (c) Rmax = 8 and R m in = 4, the imageries are subjected to 


proposed hybrid fractal image compression (HFIC) 
scheme. The Compression Ratio (CR) and Peak Signal to 
Noise Ratio (PSNR) values for the both color and gray 
images of Standard Lena, Medical MRI , pepper ,god hills 
& Zelda are determined and obtained results are tabulated 
in Table 1 and 2 respectively. The uniqueness of the 
proposed method is, it gives good comparable results 
when applied on both color and grey imageries 
The computational cost of the proposed method on Intel 
core i3 processor with 2 GB RAM is about 90 seconds for 
standard color image of Lena with size of 
5 12x5 12. Compression Time for the God hills image with 
size 600x800 is 101 seconds. When applied to compress 
the Medical images with size 512x512 time taken is about 
79 seconds. For compressing Zelda images with size 
512x512 the proposed algorithm takes about 79 seconds. 
It is observed from the table 1 that out of the three variable 
range methods, size R ma x = 16 and R m i n =8 show better 
performance in CR and R ma x = 8 and R m in=4 show better 
PSNR for test colour images. Further Table 2 indicates 
comparable performance in the CR and PSNR values for 
test grey images. 


Table. 1: CR and PSNR of Test Color Images 


s 

N 

O 

Test 

Images 

Rma X =16 and 

Rmin= 4 

Rmax =16 and 

Rmin = 8 

Rmax 8 and Rmi 

=4 

C.R 

PS 

NR 

C.R 

PSNR 

C.R 

PSN 

R 

1 

Lena 

Image 

3.09 

27.63 

17.85 

24.95 

3.27 

32.5 

6 

2 

God 

Hills 

Image 

3.96 

19.75 

16.75 

20.57 

3.18 

31.2 

5 

3 

Pepper 

Image 

3.84 

20.75 

18.95 

21.26 

3.28 

30.2 

5 

4 

Medica 
1 Image 

3.56 

21.25 

17.38 

22.49 

3.09 

31.2 

3 

5 

Zelda 

image 

3.78 

19.57 

18.83 

23.54 

3.12 

29.9 

6 


Table.2: CR and PSNR of Test Gray Images 


S 

n 

0 

Test 

Images 

Rmax=16 and 

Rmin=4 

Rmax =16 and 
Rmin=8 

Rmax 8 and Rmi 

= 4 

C.R 

PS 

NR 

C.R 

PS NR 

C.R 

PSNR 

1 

Lena 

Image 

3.01 

28.85 

17.54 

25.58 

3.04 

31.15 

2 

God 

Hills 

Image 

3.19 

22.54 

18.75 

24.36 

3.37 

31.29 

3 

Pepper 

Image 

3.35 

21.52 

19.75 

23.49 

3.59 

30.26 

4 

Medica 
1 Image 

3.89 

23.16 

17.23 

22.51 

3.42 

32.81 

5 

Zelda 

image 

3.64 

20.19 

18.75 

23.49 

3.19 

31.56 


The original and reconstructed color test images for the 
variable range block size R ma x = 16 and R m in = 8 are shown 
in Figure.4 to 8. It may be seen from the figures that the 
using variable range size for all different category images 
significantly good quality of reconstructed images is 
maintained .Similarly the results of the grey level images 
are shown in Figure. 9 to 13. 
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(b) 


(b) 



Figure .4: Results of Lena color image a) Lena Image b) re-constructed 
Lena image 



Figure 6: Results of Pepper color image a) original Pepper image b) re- 
constructed Pepper image 


(a) 


(a) 




(b) 

Figure 5: Results of God Hills color image a) original God Hills b) re- 
constructed God Hills image 


(b) 
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Figure. 7: Results of Medical MRI color image a) original Medical 
MRI image b) re- constructed Medical MRI image 







Graphics, Vision and Image Processing Journal, ISSN 1687-398X, Volume 16, Issue 2, ICGST LLC, Delaware, USA, Dec. 2016 



(b) 


Figure. 8: Results of Zelda color image a) original Zelda b) re- 
constructed zelda image 





(a) 





(b) 


(b) 


Figure. 9: Results of Lena gray image a) original lena gray Image b) re 
constructed lena image 


Fig. 1 1 : Results of pepper gray image a) original gray pepper Image b) 
re constructed pepper image 
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5.1 Comparison with other Existing Methods 

The proposed Hybrid Fractal Image Compression (HFIC) 
method obtained values for CR &PSNR are compared 
with the other existing methods like standard algorithm of 
quad tree technique[32], fractal with quad tree & Discrete 
Cosine Transform technique[33], and Fractal image 
compressions based on variable size range block[34] in 
Table 3 ,4 and Table 5,6 for color & grey images. 

The CR and PSNR values obtained from the proposed 
method in both color and gray level image for variable 
range block size R ma x = 16and R m in = 8 are comparatively 
better as depicted in figures 14,15, 16. 


Table 3: CR of the proposed method and other available methods of 
color images 


O 2 

Test Images 

Standard algorithm of 
Quadtree Technique 
(SAQT) 

Fractal with Quadtree 
and DST 

T echnique(FQ&D ST) 

FIC based on variable 
size range block 
method (FIC VRB) 

Proposed HFIC 

1 

Lena Image 

10.35 

14.03 

13.9 

17.85 

2 

Pepper 

Image 

9.75 

13.75 

14.5 

15.63 

3 

Zelda 

Image 

9.63 

14.06 

15.75 

15.48 

4 

God Hills 
Image 

9.98 

13.86 

15.25 

16.75 

5 

MRI Image 

10.25 

13.94 

14.96 

17.53 


Table 4: CR of the proposed method and other available methods of 
gray level images 


S,No 

Test Images 

Standard Algorithm of 
Quadtree 

Technique(SAQT) 

Fractal with Quadtree 
and DST 

Technique(F Q&DST) 

FIC based on variable 
size range block 
method(FIC VRB) 

Proposed 

HFIC 

1 

Lena Image 

9.72 

12.32 

13.6 

17.54 

2 

Pepper 

Image 

8.86 

12.45 

16.52 

18.75 

3 

Zelda 

Image 

9.1 

13.24 

17.32 

19.23 

4 

God Hills 
Image 

8.96 

12.68 

16.46 

19.64 

5 

MRI Image 

9.12 

12.73 

17.54 

18.94 


(b) 

Fig. 13: Results of Zelda gray image a) original gray zedla Image b) re 
constructed Zelda image 
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Table 5: PSNR of the proposed method and other available methods of 
color images 


S,No 

Test Images 

Standard 
Algorithm of 
Quadtree 

Tprlmimip fSAOT'l 

Fractal with 
Quadtree and DST 
Technique 

FIC based on 
variable size range 
block method(FIC 

Proposed 

HFIC 

1 

Lena Image 

25.64 

30.83 

25.9 

32.56 

2 

Pepper 

Image 

23.46 

29.63 

19.75 

32.56 

3 

Zelda 

Image 

22.15 

30.25 

18.64 

33.54 

4 

God Hills 
Image 

24.15 

31.64 

16.75 

34.56 

5 

MRI Image 

24.65 

30.49 

17.46 

32.49 


Table 6: PSNR of the proposed method and other available methods of 
gray level images 


S,No 

Test Images 

Standard Algorithm of 
Quadtree 

Technique(SAQT) 

Fractal with Quadtree 
and DST Technique 
(FQ&DST) 

FIC based on variable 
size range block 
method(FIC VRB) 

Proposed 

HFIC 

1 

Lena Image 

24.86 

30.18 

26.5 

31.15 

2 

Pepper 

Image 

2.75 

30.15 

17.65 

31.56 

3 

Zelda 

Image 

23.65 

29.45 

18.56 

33.57 

4 

God Hills 
Image 

24.51 

27.56 

19.43 

34.56 

5 

MRI Image 

23.49 

28.64 

19.54 

33.72 




Figure 15: Comparison chart of CR of the proposed method with other 
existing methods when applied on grey test images 


r 


>v 



Figure 16: Comparison chart of PSNR of the proposed method 
with other existing methods when applied on color test images 


6. Conclusions & Future scope 

In this paper, experimental work is done to generate hybrid 
fractal image compression using variable range block size. 
The experimental results demonstrate that the variable 
range block size of R ma x = 16 and R m in =8 gives better results 
compare to range blocks. It is clear that the proposed 
method obtains best image quality with low 
computational cost for both color and grey level images 
giving approximately same CR and PSNR values. When 
compare with the fixed block size method of fractal 
compression schemes, the present method exhibits higher 
compression ratio and PSNR values for all types of 
imageries. The proposed method is tested for different 
categories of images like standard color images, natural 
image, God Hill images, pepper image, MRI images and 
Zelda images gives comparatively good results. The PSNR 
values obtained in this approach is always greater than 30. 
The average computational cost of the proposed method 
for all categories of images is about 123 seconds. 


Figure 14: Comparison chart of CR of the proposed method with other 
existing methods when applied on color test images 


Further the proposed method image metrics can be 
improved when combined with certain lossless 
compression techniques where in machine learning 
algorithms can be used for further improving 
computational time 
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